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ABSTRACT 

Comparison of quasar absorption line spectra with laboratory spectra provides a pre- 
cise probe for variability of the fine structure constant, a, over cosmological time-scales. 
We constrain variation in a in 21 Keck/HIRES Siiv absorption systems using the al- 
kali doublet (AD) method in which changes in a are related to changes in the doublet 
spacing. The precision obtained with the AD method has been increased by a factor 
of 3: Aa/a = (—0.5 ± 1.3) X 10~ 5 . We also analyse potential systematic errors in this 
result. Finally, we compare the AD method with the many- mult iplet method which 
has achieved an order of magnitude greater precision and we discuss the future of the 
AD method. 

Key words: atomic data - line: profiles - techniques: spectroscopic - quasars: ab- 
sorption lines - ultraviolet: general 



1 INTRODUCTION 

An experimental search for variation in the fine structure 
constant, a = e 2 /Tlc, is strongly motivated by many modern 
theories. In particular, theories attempting to unify gravi- 
tation with the other fundamental interactions, such as su- 
perstring theory, often require the existence of extra, 'com- 
pactified', spatial dimensions. The scale size of the extra 
dimensions is related to the values of the fundamental con- 
stants in our observable 4-dimensional subspace. Thus, evo- 
lution of these scale sizes with cosmology leads to time vari- 
ation of the fundamental constants, such as a (e.g. Forgacs 
& Horvath 1979; Marciano 1984; Barrow 1987; Damour & 
Polyakov 1994). Currently, it seems that there is no theoret- 
ical mechanism for keeping the extra dimensions constant in 
size (Li & Gott 1998) and so such theories naturally predict 
variation of a. 

Quasar (QSO) absorption line studies provide a pow- 
erful probe of such variation. Savedoff (1956) first analysed 
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doublet separations seen in galaxy emission spectra to ob- 
tain constraints on variation of a. Absorption lines in in- 
tervening clouds along the line of sight to QSOs are sub- 
stantially narrower than intrinsic emission lines and there- 
fore provide a more precise probe of a. Bahcall, Sargent & 
Schmidt (1967) first used alkali doublet (AD) spacings of 
gas seen in absorption which seemed to be intrinsic to the 
QSO. They obtained the constraint Aa/a = (a z — ao)/ao = 
(—2 ± 5) x 10 -2 at a redshift z w 1.95. Here, a z and cvo are 
the values of a at the absorption redshift, z, and in the lab- 
oratory respectively. Since then, several authors (e.g. Cowie 
& Songaila 1995; Varshalovich, Panchuk & Ivanchik 1996) 
have applied the AD method to doublets of several species 
(e.g. C IV, Sin, Si IV, Mgll, Aim etc.) arising from interven- 
ing absorption clouds at significantly lower redshift than the 
background QSO. 

The most recent and stringent constraint using the AD 
method was obtained by Varshalovich, Potekhin & Ivanchik 
(2000, hereafter VPI00) using the Si iv A1393/1402 doublet. 
A change in a will lead to a change in the doublet separation 
given by (VPI00, correcting a typographical error) 



Aa/a — — 
' 2 



(AA/A), 
(AA/A)„ 



(1) 



Here, (AA/A) Z and (AA/A)o are the relative doublet sepa- 
rations in the absorption cloud (at redshift 2) and in the 
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laboratory respectively and c r ~ 1 is a constant taking into 
account higher order relativistic corrections. From 16 ab- 
sorption systems (towards 6 QSOs) they obtained a mean 
Aa/a = (-4.6 ± 4.3) x 10~ 5 using a line fitting method. 
The error quoted here is statistical only. VPI00 augment this 
with an additional systematic error term, ±1.4x 10 -5 , due to 
uncertainties in the laboratory doublet separation, (AA/A)o, 
assumed in their analysis (see Table |lj) . This estimate of the 
potential systematic error seems optimistic considering that 
the error in the laboratory wavelength separation was pre- 
viously quoted at S(AX/X)o ~ 1mA (Ivanchik, Potekhin & 
Varshalovich 1999). From equation [l], the corresponding sys- 
tematic error in Aa/a is 



\5{Aa/a) 



c r S(A\/\) 
' 2 (AA/A) 



5 x 10" 



(2) 



which is more consistent with the value reported by Ivanchik 
et al. (1999) (~8x 1CT 5 ) than by VPI00. 

There are three ways in which the above constraints can 
be significantly improved without the need for a much larger 
sample: 

(i) Improved spectral resolution: Many of the spectra 
used by VPI00 have FWHM ~ 20kms _1 (e.g. Petitjean, 
Rauch & Carswell 1994; Varshalovich et al. 1996) while 
many absorption systems are known to contain Siiv lines 
with 6-parameters (i.e. Doppler widths) ~ Skms" 1 (e.g. 
Outram, Chaffee & Carswell 1999). FWHM ~ Tkms" 1 is 
now routinely used for QSO observations with, for example, 
the Keck/HIRES and the VLT/UVES. Therefore, such spec- 
tra offer the opportunity to significantly improve constraints 
on Aa/a. 

(ii) Improved signal-to-noise ratio (SNR): Many of the 
spectra used by VPI00 have SNR ~ 15 per pixel. With the 
advent of 8-10 m telescopes, significantly greater SNR can 
be achieved (e.g. Prochaska & Wolfe 1996, 1997, 1999). 

(iii) Improved laboratory wavelength measurements: 
Griesmann & Kling (2000) have increased the absolute pre- 
cision of the laboratory wavelengths of Siiv A1393 and 1402 
by more than two orders of magnitude. We compare the 
new values with those used by VPI00 (Kelly 1987; Morton 
1991, 1992) in Table 1. Note that even our conservative es- 
timate of the systematic error on Aa/a in equation ^| is 
too small. The new measurements imply that the VPI00 re- 
sult in equation |l| should be corrected by ~ 11.3 x 10 -5 
to Aa/a w (7 ± 4) x 10~ 5 . If the error in the differ- 
ence between the doublet wavelengths is equal to the ab- 
solute uncertainty (i.e. 4 x 10 -5 A) then a precision limit of 
8 (Aa/a) ~ 0.2 x 10 -5 can now be reached. This should be 
true of Griesmann & Kling's measurement since statistical 
errors in the line positions dominated their error budget (Ulf 
Griesmann, private communication). 

We apply the above improvements in the present work 
with the aim of increasing both the accuracy and precision 
of Aa/a as measured using the Siiv AD method. The paper 
is organised as follows. Section ^ describes our observations, 
analysis methods and results. We investigate potential in- 
strumental and astrophysical systematic effects in Section 
|i| We summarise our work in Section ^ and compare our 
results with others in the literature. We also compare the 
AD method with another method, the many-multiplet (MM) 



method (Dzuba, Flambaum & Webb 1999a, b; Webb et al. 
1999), which offers an order of magnitude greater precision. 



2 DATA, ANALYSIS AND RESULTS 
2.1 Keck/HIRES observations 

All our QSO spectra were obtained at the Keck I 10 m tele- 
scope on Mauna Kea with the HIRES facility (Vogt et al. 
1994) over several observing runs from 1994 to 1997. The 
QSOs were generally quite faint (mv ^ 19.0) so several ~1— 
2 hour exposures were co-added for each object. Most of the 
data were reduced using the HIRES data reduction pack- 
age written by T. Barlow, makee. This package converts 
the two-dimensional echelle images to fully reduced, one- 
dimensional, wavelength-calibrated spectra. 

Thorium-Argon (ThAr) calibration spectra were taken 
before and after the QSO exposures and co-added to pro- 
vide a calibration spectrum. ThAr lines were selected and 
centroided to provide a wavelength solution. Some of the 
spectra were reduced when makee had no wavelength cal- 
ibration facility. In these cases, wavelength calibration was 
carried out using irafQ routines. Spectra not reduced in 
makee were fully reduced within iraf. la error arrays were 
generated assuming Poisson counting statistics. We fit con- 
tinua to regions of each spectrum containing either or both 
of the Siiv doublet transitions by fitting Legendre polyno- 
mials to ~ 500 kms" 1 sections. Full details of the reduction 
procedures can be found in Prochaska & Wolfe (1996, 1997, 
1999). Outram et al. (1999) have also kindly contributed 
their spectrum of Q 1759+75 taken in July 1997. 

Our sample comprises 21 Siiv absorption systems (to- 
wards 8 QSOs) over a redshift range z = 2-3 (mean z = 2.6). 
The SNR per pixel ranges from 15-40 with most spectra hav- 
ing SNR ~ 30 and FWHM <, 7.5 kms" 1 (R = 34000). We 
provide an example absorption system in Fig. 1. 



2.2 Analysis 

Although equation [l] is a simple approach to the specific 
case of an alkali doublet, a more general approach is to 
write down the energy equation for an individual transi- 
tion, within any multiplet and for any species. Dzuba et al. 
(1999a,b) and Webb et al. (1999) suggested the convenient 
formulation 



loq + qix + q2y , 



(3) 



where ui z is the wavenumber in the rest-frame of the cloud, 
at redshift z, in which a z /ao may not equal unity, cjq is 
the wavenumber as measured in the laboratory and x and 
y contain the information about Aa/a: x = (^j) — 1, y = 

(fp) 4 — 1- The qi and q2 coefficients represent the relativistic 
corrections to the energy for a particular transition. The 
qi coefficients are typically an order of magnitude larger 
than the q-i coefficients and so it is the relative magnitudes 



t iraf is distributed by the National Optical Astronomy Obser- 
vatories, which are operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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Table 1. Atomic data for the SilvA1393 and 1402 lines. We give the ground and excited state configurations and compare the new 
laboratory wavelength measurements with those used by VPI00. We also give the oscillator strengths, /, used in our profile fitting 
algorithm and the relativistic coefficients qi and q 2 used in equation H. 



Transition 


Ground 


Upper state 


Old A /A a 


New A /A 6 


New u>o 1 cm 1 b 


f a 9i/cm lc Q2 /cm lc 


SiivA1393 
Si iv A1402 


2p 6 3s 2 S 1/2 


2p 6 3p 2 P 3/2 
2p 6 3p 2 P 1/2 


1393.755(6) 
1402.770(6) 


1393.76018(4) 
1402.77291(4) 


71748.355(2) 
71287.376(2) 


0.5140 766 48 
0.2553 362 -8 



"Martin & Zalubas (1983) and Kelly (1987); 6 Griesmann & Kling (2000); c Dzuba, Flambaum & Webb (1999b) 





Velocity (km/s) 

Figure 1. Si IV absorption system at z = 2.530 towards the QSO 
2348-14. The data have been normalized by a fit to the continuum 
and plotted as a histogram. Our Voigt profile fit (solid curve) 
and the residuals (i.e. [data] — [fit]), normalized to the \a errors 
(horizontal solid lines), are also shown. The tick-marks above the 
continuum indicate individual velocity components. 



of q\ for different transitions that characterizes our ability 
to constrain Aa/a. This equ ation forms the basis of the 
MM method (see Section 4.1). In the case of a single alkali 
doublet and Aa/a <C 1, equation ^ reduces to equation [l] 
with 



<5gl + 5q 2 
5qi + 2Sq 2 



(4) 



where Sqi and 8q 2 are the differences between the qi and q 2 
coefficients for the doublet transitions. The values for q\ and 
q 2 for Si IV A1393 and 1402 have been calculated by Dzuba et 
al. (1999b) and are shown in Table 0. For the Siiv doublet, 
c r w 0.9 as used by VPI00. 

Our technique is based on a simultaneous \ 2 minimiza- 
tion analysis of multiple component Voigt profile fits to the 
Si IV absorption features in the QSO spectra. Consider a 
QSO spectrum containing a single velocity component of 
a specific transition. Three parameters describe the Voigt 
profile fit to such a component: the column density N, the 
Doppler width or fa-parameter and the redshift z of the ab- 
sorbing gas cloud. For the present case, we add another free 
parameter to the fit: Aa/a. 



We have used the program vpfit^ (Webb 1987) to fit 
absorption profiles to the spectra. We have modified vpfit 
to include Aa/a as a free parameter. Parameter errors can 
be calculated from the diagonal terms of the final parameter 
covariance matrix (Fisher 1958). The reliability of these er- 
rors has been confirmed using Monte Carlo simulations of a 
variety of different combinations of transitions and velocity 
structures. 

vpfit imposes a cut-off point in parameter space such 
that very weak velocity components are removed from the 
fit when they no longer significantly affect the value of \ 2 - 
Conceivably, dropping even very weak components could af- 
fect our determination of Aa/a. Therefore, in such cases 
we observed the trend in the values of (Aa/a); at each it- 
eration i of the minimization routine to see if this trend 
was significantly altered due to line dropping. If components 
were dropped during a fit then we also re-ran the vpfit al- 
gorithm, keeping the dropped components by fixing their 
column density at the value just before they were dropped 
from the original fit. No cases were found where the values 
of Aa/a from the different runs differed significantly. 

We impose several consistency checks before we accept 
a value of Aa/a. Firstly, the value of \ 2 P er degree of free- 
dom must be ~1. Secondly, we calculated Aa/a for each 
absorption system using a range of different first-guess val- 
ues for Aa/a to ensure that the vpfit algorithm was finding 
global minima in the \ 2 parameter space. 



2.3 Results 

We present our results in Table ^. We show our values of 
Aa/a for each absorption cloud together with the la er- 
ror bars derived from the vpfit algorithm. To illustrate the 
distribution of Aa/a over cosmological time, we plot these 
results in Fig. 2 as a function of fractional look-back time to 
the cloud and the absorption cloud redshift, z. The weighted 
mean of the sample is Aa/a = (—0.5 ± 1.3) x 10~ 5 . Our re- 
sults show a 3.3-fold increase in precision over the VPI00 
result. The reduced \ 2 about the weighted mean value is 
0.95, which gives no evidence to suggest that we may have 
incorrectly estimated the individual la errors on Aa/a. We 
checked that the distribution of Aa/a values shows no gross 
deviation from a Gaussian distribution (although with only 
21 data points this is not a rigorous check). We also note 
that the unweighted (Aa/a = (-0.12 ± 1.4) x 10" 5 ) and 



See 



http://www.ast.cam.ac.uk/~rfc/vpfit.html for details 



about obtaining VPFIT. 
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Table 2. The raw results from the \ minimization procedure. 
For each QSO sight line we identify the QSO emission redshift, 
2 cm , the absorption cloud redshift, 2 abs and the value of Aa/a 
for each absorption cloud. 



Redshift 



Object 



Aa/a/10- 



0100+13 


2.68 


2.299 


-3.05 + 7.30 






2.309 


8.96 ± 10.39 


0149+33 


2.43 


2.065 


-4.47 + 16.81 






2.140 


-12.93 ± 12.69 






2.204 


-5.89 + 6.68 


0201+36 


2.49 


2.457 


-4.50 + 3.57 


0347-38 


3.24 


2.810 


3.04 ± 7.43 






2.899 


-10.90 ± 10.87 






3.025 


-4.05 + 6.78 


1759+75 


3.05 


2.624 a 


-5.98 + 2.98 






2.849 a 


2.56 + 3.29 






2.849 


2.74 + 3.88 






2.911 


6.00 ± 9.89 






2.911 a 


0.65 + 4.75 


2206-20 


2.56 


2.014 


-3.84 + 6.59 






2.128 


-3.48 ± 7.41 


2231-00 


3.02 


2.641 


0.20 ± 10.65 






2.986 


4.23 ± 17.74 


2348-14 


2.94 


2.279 


10.38 + 6.05 






2.530 


3.76 + 6.53 






2.775 


14.04 + 7.10 



"These absorbers contributed by Outram, Chaffee 
(1999) 



Carswcll 



weighted means are consistent with each other, again sug- 
gesting no substantial deviation from Gaussianity. 



3 POSSIBLE SYSTEMATIC ERRORS 

Although our result is statistically insignificant, the preci- 
sion obtained is approaching the level at which systematic 
errors begin to become important. We have already dis- 
cussed how the VPI00 result was severely effected by er- 
rors in the laboratory wavelengths of SilvA1393 and 1402. 
We must therefore examine possible systematic errors in our 
results. We have considered a range of instrumental and as- 
trophysical systematic errors particular to the MM method 
in Murphy et al. (2001b, hereafter MOlb). Below we discuss 
some of these in the context of the Si IV AD method. 



3.1 Wavelength mis-calibration 

The wavelength calibration of the QSO CCD images is done 
by comparison with Thorium-Argon (ThAr) lamp exposures 
taken before and after the QSO frames. We consider errors 
in the laboratory wavelengths of the selected ThAr lines to 
be negligible as their relative values are known to a similar 
precision as our values of uio (Palmer & Engleman Jr. 1983) . 
However, ThAr line mis-identifications may lead to serious 
mis-calibration of the wavelength scale over large wavelength 
regions. Such mis-identifications can also be applied to the 
rest of the spectra taken over the same series of observations 
(i.e. applied to the spectra of many QSOs). Thus, a priori, 
we cannot rule out the possibility that the wavelength scale 
has been set improperly in this process and that a systematic 



20 



H =68 km/s/Mpc 
- (0 M ,Q A ) = (0.3,0.7) 
t = 13.86 Gyr 



a 



-20 



0.76 



0.78 0.8 0.82 0.84 
Fractional look-back time 



Figure 2. Our values of Aa/a for each absorption cloud plotted 
against fractional look-back time (using the specified cosmological 
model) and redshift. The weighted mean of the sample is Aa/a = 
(-0.5 + 1.3) X 10" 5 . 



shift in a has not been mimicked. Such a potential effect 
needs careful investigation. 

We designed a direct test for wavelength mis-calibration 
in MOlb: we treated the ThAr emission lines in the calibra- 
tion spectra with the same analysis used to derive Aa/a 
from the QSO absorption lines in the object spectra. We 
have applied this method to the ThAr spectra corresponding 
to the absorption systems in our sample. For each absorption 
system we derive values of (Aa/cv)ThAr from corresponding 
sections of ThAr s pect ra using the same analysis technique 
outlined in Section 2.2 with the following modifications: 



(i) We fit Gaussian profiles to the ThAr emission lines 
(instead of Voigt profiles which were fitted to the QSO ab- 
sorption lines). 

(ii) We use the literature wavenumbers of the selected 
ThAr lines (Palmer & Engleman Jr. 1983) for loq in equation 

However, we use the qi and q2 coefficients of the Si IV lines. 
For example, if we select a ThAr line lying in the same region 
of the spectrum as the SilvA1393 line of interest, then we 
assign the qi and q2 coefficients for Siiv A1393 to that ThAr 
line. 

(iii) In MOlb we noted that the x 2 f° r a fit to a set of 
ThAr lines was 3> 1 (in the present case, this 'set' contains 
two ThAr lines lying in the same spectral regions as the cor- 
responding Si IV lines). This was attributed to weak blends 
in the selected ThAr lines. This was also observed in the 
present analysis. In this case, we cannot use the lcr error on 
(Aa/a)ThAr generated by vpfit. Therefore, we select sev- 
eral different sets of ThAr lines and find (Aa/a)ThAr for 
each set. The error in (Aa/a)-rhAr is then taken as the error 
in the mean value for the different sets. 

Our results are shown in Fig. 3 where we compare them 
with the results from the QSO data. It is clear that the QSO 
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Figure 3. Comparison of QSO and ThAr results. The top panel 
shows the QSO results and the middle panel shows the ThAr data 
on the same scale. The weighted mean is (-6 ± 9) x lCT 7 . The 
lower panel is an expanded view of the ThAr results. Note that 
some ThAr data was not available. Clearly, this does not affect 
our conclusions. 



results have not been significantly affected in any systematic 
way. The weighted mean value of (Aa/a)xhAr = (— 6±9) x 
10 -7 . The scatter in the data is larger than what we expect 
on the basis of the individual error bars. This could be due to 
either real (low level) mis-calibrations or to the line blending 
discussed above. 

The results in Fig. 3 are also relevant for possible vari- 
ations in instrumental profile (IP) asymmetry along and 
across echelle orders. Valenti, Butler & Marcy (1995) have 
determined the HIRES IP for several positions along a single 
order and they find that the IP is indeed asymmetric and 
that the asymmetry varies slightly along the echelle order. 
However, Fig. 3 shows that any asymmetry variations have 
not systematically affected our determination of Aa/a. 



3.2 Systematic line blending with unknown 
species 

There may have been weak, interloping, unresolved lines 
which, if the interloping species were in the same absorption 
cloud, could have produced a shift in the fitted line wave- 
lengths of all velocity components of one or both Si IV tran- 
sitions. We distinguish between random blends and system- 
atic blends. Random blends may occur if many absorption 
clouds at different redshifts intersect the line of sight to a 
single QSO. A systematic blend will occur when two species 
are in the same cloud and have absorption lines with similar 
rest-wavelengths. Such an effect could mimic a systematic 
shift in a. 

We have searched atomic line databases (Moore 1971; 
the Vienna Atomic Line Database (VALD) - Piskunov et al. 
1995 and Kupka et al. 1999) for transitions of any species 



which may blend with SilvA1393 or 1402. The search was 
restricted to transitions from the ground state with rest 
wavelengths A such that | A — Ao | < 0.2 A. This is a con- 
servative upper limit adopted from simulations of typical 
blends in MOlb. We did not identify any potential blends 
with SilvA1402 satisfying this criteria. For SilvA1393 we 
found one potential interloper: BelA1393 at A = 1393.804 A 
(VALD). This transition has an oscillator strength of / = 



2.080 x 10" 



The Be I A2349 transition is nearly three or- 



ders of magnitude stronger than this (/ = 1.698) and, to 
our knowledge, has never been detected in QSO absorption 
spectra. This safely rules out Be I A1393 as a candidate in- 
terloper. 

In summary, we have found no known atomic transi- 
tions which may blend with either of the Siiv doublet lines. 
We have not considered molecular transitions as potential 
interlopers but consider this possibility to be unlikely (see 
MOlb for discussion). 



3.3 Differential isotopic saturation 

Si has three naturally occurring isotopes, 28 Si, 29 Si and 30 Si, 
with terrestrial abundances in the ratio 92.23 : 4.68 : 3.09 
(Rosman & Taylor 1998). Thus, each absorption line will 
be a composite of absorption lines from all three isotopes. 
We are not aware of any experimental determinations of the 
spectral separations between the isotopic components. Thus, 
the values of luq in Table |l| (Griesmann & Kling 2000) are 
composite wavenumbers only. These wavenumbers will only 
strictly be applicable in the optically thin regime (linear part 
of the curve of growth). As the column density increases, the 
strongest isotopic component begins to saturate and the line 
centroi d wi ll shift according to the natural abundances (cf. 
Section 



3.4) of the other isotopes. 



Estimates of the isotopic separations for the Siiv dou- 
blet have been made by W. R. Johnson (private communica- 
tion) and these are summarised in Table [| From these sep- 
arations we have estimated the absolute isotopic wavenum- 
bers (coo in Table |^) such that the abundance- weighted mean 
of the isotopic wavenumbers equals the composite values in 
Table |l| The error in these results is Ss 10% and comes from 
errors in the isotopic separations. 

If we consider a Si IV doublet with a single velocity com- 
ponent then, as the A1393 line begins to saturate, the line 
centroid will shift bluewards. The effective doublet separa- 
tion will continue to increase but this increase will eventually 
slow and reverse as the A1402 transition begins to saturate. 
It is clear that this is a potentially important systematic 
effect for the Siiv doublet. Therefore, we have used the val- 
ues of uo in Table ^| to re-calculate Aa/a for our sample 
of Si IV absorbers. Only 5 systems clearly contain saturated 
components and Aa/a in these systems displays the ex- 
pected behaviour: Aa/a becomes more negative when the 
isotopic structures are included. Overall, the weighted mean 
for our sample becomes Aa/a = (—0.8 ± 1.3) x 10 -5 . 



3.4 Isotopic abundance variation 

Timmes & Clayton (1996) suggest that the abundance of 
29 Si and 30 Si should decrease with decreasing metallicity. In 
general, damped Lyman-a systems (DLAs), which some of 
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Table 3. The isotopic structures of SilvA1393 and 1402 calculated by W. R. Johnson (private communication). We show the total 
isotopic separation, Auj^ ot , and the contribution from the specific (Aojg Pec ) and volume (AoJq 01 ) isotopic shifts. The final wavenumber, 
luq, is such that the abundance- weighted mean equals the corresponding composite value (see Table [l]) . The last column shows the relative 
abundance of each isotope (Rosman & Taylor 1998). 



Transition 


m/amu 


A^/cm" 1 


Au^/cm- 1 


Aa# ol /10- 3 cm" 1 


UJ0 1 cm 1 


A /A 


% 


SilvA1393 


28 








71748.344 


1393.76039 


92.23 




29 


0.1044 


0.0585 




71748.478 


1393.75779 


4.68 




30 


0.2017 


0.1130 


4.413 


71748.545 


1393.75649 


3.09 


SiivA1402 


28 








71287.365 


1402.77313 


92.23 




29 


0.1042 


0.0583 




71287.499 


1402.77049 


4.68 




30 


0.2013 


0.1126 


4.418 


71287.566 


1402.76917 


3.09 



our Si IV doublets are associated with, are thought to have 
metallicities Ss — 1 at high redshift (e.g. Prochaska & Wolfe 
2000). If the isotopic abundances of Si at redshifts z — 2- 
3 differ substantially from their terrestrial values then this 
may mimic a shift in a. 

The dominant component of the isotopic shift is the 
mass shift: Alu oc ui/m 2 for m the atomic mass. Thus, for 
a single ion, the mass shift is degenerate with the redshift 
parameter fitted to the velocity components of an absorp- 
tion system. That is, it will have no affect on our values 
of Aa/a. The specific isotopic shifts for the Siiv doublet 
have been calculated by W. R. Johnson (private communi- 
cation) and are given in Table ^| We also show the volume 
isotopic shift for 30 Silv as calculated by V. A. Dzuba (pri- 
vate communication) and these are clearly negligibly small. 
The specific shift is not degenerate with redshift but only 
at a negligible level. Thus, even if the isotopic ratios of 29 Si 
and 30 Si are much smaller in our Si IV absorbers, this will 
have a negligible affect on Aa/a. 



3.5 Hyperflne structure effects 

Hyperfine splitting can lead to saturation effects similar 



to the isotopic saturation effects discussed in Section 3.3 
The transitions of the 29 Si isotope will experience hyperfine 
splitting since the neutron number is odd. As with the iso- 
topic structure, the magnitude of hyperfine splitting has not 
been measured. However, V. A. Dzuba (private communica- 
tion) has estimated the splitting and finds it to be a factor 
of ~2 times smaller than the isotopic splitting. This may 
also be estimated by comparison with the hyperfine struc- 
ture of 25 Mg since the Si IV and Mg n doublets have similar 
ground and excited state wavefunctions and the magnetic 
moments of Mg and Si are similar. Drullinger, Wineland & 
Bergquist (1980) have measured the hyperfine splitting of 
the 25 MgnA2796 transition: the hyperfine splitting is ~2 
times smaller than the isotopic separations. It is therefore 
clear that saturation effects will be negligible. 



3.6 Atmospheric dispersion effects 

As discussed in MOla, 5 of our absorption systems were ob- 
served without the use of an image rotator (i.e. they were 
observed before August 1996 when a rotator was installed 
on HIRES). The QSO light is dispersed across the spectro- 
graph slit by the atmosphere and, if the slit is not maintained 
perpendicular to the horizon by an image rotator, light of 



different wavelengths will enter the spectrograph at slightly 
different angles. In MOlb we discussed two effects resulting 
from this: (i) a bulk stretching of the spectra relative to the 
ThAr calibration frames and (ii) a distortion of the PSF due 
to truncation of the dispersed seeing disc on either side of the 
slit jaws. We do not consider (ii) any further since the effect 
on Aa/a found in MOlb was small and since the distortion 
should be very similar at both Si IV doublet wavelengths. We 
consider effect (i) below. 

The optical design of HIRES at the time at which our 
affected data was taken was such that 9 = £ for 6 the angle 
of the slit to the vertical and £ the zenith angle of the object 
(Tom Bida, Steve Vogt, personal communications). Consider 
two wavelengths, Ai and A2 (A2 > Ai), falling across the 
spectrograph slit. If we were to measure the spectral sepa- 
ration between these two wavelengths on the CCD, AA', we 
would find that 



AA' « A 2 - Ai + 



aAij) sin 9 



(5) 



where a is the CCD pixel size in Angstroms, S is the 
projected slit width in arc seconds per pixel (for HIRES, 
S — 0'.'287 per pixel) and Aip is the angular separation (in 
arc seconds) of Ai and A2 at the slit. Aip is a function of the 
atmospheric conditions along the line of sight to the object 
and can be approximated using the refractive index of air 
at the observer and the zenith distance of the object (e.g. 
Filippenko 1982). Note that equation ^ is only valid if the 
seeing and tracking error are zero. These two effects will 
reduce the measured AA'. 

For the Siiv doublet, a typical value for Aip is only w 
0'.'007 (using atmospheric conditions typical of Mauna Kea, 
an absorption redshift z — 2.6 and ip = 30°, which is typical 
of our sample, see MOlb). This implies AA' — (A2 — Ai) » 
3.7 x 10 -4 A with a — 0.03 A per pixel and will mimic a shift 
in a, Aa/a ~ +0.66 x 10 -5 . The sign of this correction is 
such that removing the stretching implied by equation ^ 
from our data will result in a more negative Aa/a. The 
correction is substantial for an individual system so we have 
applied it to the 5 affected systems to obtain an upper limit 
on the overall affect on Aa/a. Our weighted mean value of 
Aa/a for the entire sample becomes (—0.6 ± 1.3) x 10 . 



4 DISCUSSION 
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4.1 Comparison with other results 

We have increased the precision of Aa/a as measured with 
the AD method by at least a factor of 3.3 over that of 
VPIOO. The increase in precision (and accuracy) is due to 
improved data quality (higher spectral resolution and SNR) 
and greatly improved laboratory wavelength measurements 
(Griesmann & Kling 2000). 

Despite this increase in precision, the MM method can 
provide a further increase of up to an order of magnitude 
with similar quality data. The MM method was first pro- 
posed by Dzuba, Flambaum & Webb (1999a,b) and Webb 
et al. (1999) (see also Dzuba et al. 2001). It is based on equa- 
tion tfl and the use of many transitions, from many different 
multiplets in different ions. The qi coefficients of transitions 
in heavy (m ~ 50) ions are typically an order of magnitude 
larger than those for light ions (m ~ 25). Depending on the 
nature of the ground and excited state wavefunctions, the 
qi coefficients for some transitions may also be of opposite 
sign, even for transitions of the same ion. Therefore, the MM 
method has the following advantages over the AD method: 

• The MM method allows us to probe an effect ~10 times 
larger since we can compare total relativistic corrections (not 
just the spin orbit coupling probed by the AD method). 

• In principle, all transitions appearing in the QSO spec- 
tra may be used and so the number of constraints on Aa/a 
is increased. 

• In practice, having a larger number of transitions leads 
to much better constraints on the velocity structure of the 
absorbing clouds. 

• Comparison of transitions with positive and negative qi 
coefficients minimizes the systematic effects in Section |j| 

Webb et al. (1999) demonstrated the MM method and 
found tentative evidence for a smaller a at redshifts z ~ 1: 
Aa/a = (-1.09 ± -0.36) x 10 -5 . We have extended that 
work in MOla and Webb et al. 2001, analysing 21 high red- 
shift DLAs and re-analysing 28 lower redshift Mg/Fe sys- 
tems from the Webb et al. (1999) data set. We confirmed the 
Webb et al. result and found similar results for the DLAs in 
Murphy et al. (2001a) and Webb et al. (2001). For the en- 
tire sample (0.5 < z < 3.5), Aa/a = (-0.72 ±0.18) x 10~ 5 . 
We analysed potential systematic errors in MOlb and found 
that atmospheric dispersion effects and isotopic abundance 
variation may have led us to find a Aa/a that was too pos- 
itive (i.e. once these two effects are removed from our data, 
we find a more negative Aa/a). The results of MOla and 
MOlb indicate possible variation in a. Even if these results 
are considered to be an upper limit on Aa/a, they are by 
far the most restrictive limits at high redshift. Our present 
result of Aa/a = (— 0.5±1.3) x 10 -5 is consistent with these 
results. We review constraints on Aa/a at other redshifts 
in MOla. 



4.2 The future for the Siiv AD method 

If the AD method were to be used to reach a similar pre- 
cision as reached in MOla, a SNR J> 200 in ~50 systems 
would be required. In Section [] we found that system- 
atic errors important for the AD method can correspond 
to |Acv/a| ~ 0.3 x 10 which is of similar magnitude to 
the effect of systematic errors in the MM method (MOlb). 



Therefore, once these systematic errors are understood and 
can be reliably removed from QSO absorption data, the MM 
method is clearly the preferred method for probing Aa/a at 
high redshift. 

Despite this, a modest increase in precision could eas- 
ily be achieved with our present data if the C iv A1548/1550 
doublet laboratory wavelengths were known to high enough 
precision. Griesmann & Kling (2000) have increased the lab- 
oratory precision by an order of magnitude but the statis- 
tical uncertainty in the wavelengths is still Jj 0.04 cm -1 . If 
the laboratory precision can be increased by another order 
of magnitude, the Si IV and C IV doublet can be analysed 
simultaneously to obtain tighter constraints on Aa/a (i.e. 
the MM method could be used). 
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